Fermi National Accelerator Laboratory pT™*^~tA3°3T"" 



ON THE COMPUTATION OF CMBR ANISOTROPIES FROM 
SIMULATIONS OF TOPOLOGICAL DEFECTS 



A. Stebbins and S. Dodelson 
NASA/Fermilab Astrophysics Center, 
Fermtlab, Box 500, Batavia, IL 60510, USA 

Techniques for computing the CMBR anisotropy from simulations of topological 
defects are discussed with an eye to getting as much information from a simulation 
as possible. Here we consider the practical details of which sums and multiplica- 
tions to do and how many terms there are. 



Probably the most fruitful field of eosmological observation in the past few 
years comes from observations of the anisotropics of the Cosmic Microwave 
Background Radiation (CMBR). One of the first and the most dramatic de- 
tections was made by the DMR on the COBE satellite over 5 years ago. New 
results continue to pour in and we can expect this continue to at least until 
the Planck experiment provides us with an exquisitely detailed all-sky map in 
about a decade. Theoretical modeling of the data and their interpretation in 
terms of cosmogenic theories has generally been able to keep pace with these 
observations but we can expect the observational data to become increasingly 
more precise. There are models such as inflation where the statistical proper- 
ties of the anisotropy pattern is well described as being statistically isotropic 
Gaussian random noise where the parameters describing the Gaussian distribu- 
tion can be computed by solving linear sets of ordinary differential equations. 
For such models a complete statistical description is fairly easy to obtain and 
the predictions are routinely made made with better than 1% accuracy. In 
another class of theories, where the inhomogeneities are seeded by topological 
defects, it is not clear whether we will ever obtain this kind of accuracy. Of 
course this may be a moot point if the predictions for these models, such as 
they are, do not fit observations. In the paper we will discuss some of the 
issues involved with computing CMBR anisotropics for defect models. 

One of the fundamental differences between Gaussian models and defect 
models is that the dynamics of the former are described by linear equations 
while that the latter by nonlinear equations. Linearity combined with the as- 
sumed statistical homogeneity of the universe guarantee that one may decom- 
pose the eosmological perturbations into eigenmodes, which are eigenfunctions 
of the generators of the isometrics of the space-time, and that each eigenfunc- 
tion will evolve independently. So, for example in a flat Friedmann-Robertson- 
Walker (FRW) spacetime, the eigenmodes could be Fourier modes which are 
eigenfunctions of the generator of translation. While one might have to solve 
the linear equations for a Fourier mode numerically, one need only solve the 
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equation once for each value of |k|. For a given type of inhomogeneity, e.g. 
adiabatic or isocurvature, the only freedom in the initial condition of a given 
mode is the amplitude, and, since the equations are linear, solutions scale 
linearly with the initial amplitude. 

For topological defects the equations are non-linear so different eigenmodes 
are coupled and one does not know the evolution of one mode without taking 
into account the evolution of all the other. It is usually simpler to think 
about topological defects in real space rather than the space of eigenmodes 
(A:-space) . One reason is that equations describing the evolution of the defects 
are spatially local while they are highly non-local in /c-space. Another reason 
is that causality guarantees that the range interaction does not extend beyond 
a certain distance, the causal horizon, while the range of interaction in A;-space 
is unbounded. Causality allows one to consider a finite patch of the universe 
and ignore what goes on outside of that patch, at least for points more than a 
horizon distance away from the boundary of the patch. 

In order to determine the likelihood that the observed inhomogeneities 
in the universe could be explained by topological defects one would need a 
description of the statistical ensemble of possible topological defects configu- 
rations and their evolution. In general we cannot even find analytic solution 
to describe the evolution of specific configurations so the way one tries to de- 
termine the statistical properties is through numerical experiment. Namely 
one generates a realization of the defect configuration at some early time and 
evolves it and the matter surrounding it according to their equations of mo- 
tion. One might then try another configuration, and so on, until one has a fair 
statistical sample. It is generally believed that the statistical properties of the 
defects are ergodic so that one may also obtain a fair sample by considering a 
large enough volume of a single realization. One need only consider the obser- 
vations predicted for observers situated in different locations in the simulation 
volume. Since one needs to build up a fair sample it is important to obtain as 
many realizations of the anisotropy as possible. 



1 Geometry of Computation 



Numerous authors have used si 
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anisotropy patterns from defectsEi'D 
to examine large angle anisotropics □ 



'ects to produce realizations of 
Some of these have been used 
iJe.^others have been used to 
examine anisotropics on smaller angular scales IU'Q'Ej. The geometry used for 
these large angle numerical experiments is shown in fig One must evolve 
the defects in a very large simulation volume, as big as the present horizon 
and ideally one would like to evolve from time of recombination. Usually one 
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Figure 1: Here is shown the geometry for large angle anisotropy computations. In the inset 
the cube represents the simulation volume while the sphere gives the observer's horizon, the 
volume which contributes to the anisotropy pattern of a given observer. The shaded on the 
graph gives the range of temporal and length scales resolved by the simulation. Generally 
the most important missing scales are on the left, i.e. from early times. Some numbers 
appropriate for studying COBE scale anisotropics are given. 
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takes periodic boundary conditions for the defect configuration. However this 
periodicity does not really matter - since for most methods of laying down 
the defects one can show that the defect configuration in any given horizon 
volume is drawn from the same distribution as if one were to have simulated 
and infinite volume. The cubical box with periodic boundary condition are 
only easily implemented in the case of spatially flat FRW cosmology. Doing 
this computation for an open universe, as was done in by Pen and SpergeB, 
is actually a more difficult problem in many respects. For simplicity we will 
assume a flat FRW cosmology, which allows us to make use of the Fourier 
series representations of distributions in the simulation box, i.e. 

F(x)=^F(k)e*-. (1) 

k 

We will specifically be considering a a cubical simulation volume with N x N x 
N spatial resolution elements. 

2 Dynamic Range 

Given these general considerations let us proceed to consider specific issues 
associated with the computation of MBR anisotropics from seeds in a box of a 
fixed comoving size, L. All numerical simulations have a finite dynamical range 
in both spatial and temporal scales. There will be a finite spatial resolution 
on small scales and the size of the box itself will limit the largest accessible 
scale. Similarly there will be finite temporal resolution, usually related to 
spatial resolution by the Courant condition At < Ax/c. Here c is the speed- 
of-light which characterizes the dynamics of most defects. For this reason the 
coherence length of the defects will scale with the causal horizon, so that in 
order to properly represent the defect configuration the finite spatial resolution 
leads to a minimum starting time: ti ^ Ax. This large characteristic velocity 
of the defects also means that the defect configuration will start to "see" the 
finite size of the box when the causal horizon reaches the box size. It is desirable 
to end the simulations before this happens, i.e. tf < L/c. These considerations 
indicate that the number of temporal resolution elements is N. 

From a box at least as large as twice the horizon one can study the very 
largest angle anisotropics since, as illustrated in fig |l], one does not see the 
periodicity of the universe. The edge of ones horizon in one direction nearly 
touches the edge in the opposite direction, but one is seeing this far only at 
the earliest times when the (causal) correlation length is arbitrarily small, so 
there need be no effect of periodicity at all. The finite spatial resolution in 
such a simulation will limit the angular resolution. If most of the anisotropy is 
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Figure 2: Here is shown the geometry for a small-angle anisotropy computations. The cube 
represents the simulation volume and the cone the photons from the simulation converging 
toward the observer. The shaded region on the graph gives the scales resolved by the 
simulation. The most important missing scales are on the right and left, i.e. early times and 
late times. Some numbers appropriate for studying degree scale anisotropy are given. 
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generated at large redshifts then the angular resolution is AO ^ Ax/{L/2) ^ 
2/N, which corresponds to a largest available angular wavenumber ^max ~ N/2. 

For simulations used to compute smaller angle anisotropy, as illustrated in 
fig ^, the effects of periodicity are much more severe as one can "see" beyond 
the simulation volume. One is not really simulating a large universe as we live 
in, but rather a truly small and periodic one. In such a universe one would see 
the pattern of anisotropy to be roughly periodic on the sky. What one is hoping 
is that the statistical properties of each of these patches are similar to those of 
a similar sized patch in a non-periodic universe. The angle subtended by the 
simulation box gives you an upper limit on the angular scale you are probing, 
while the spatial resolution gives a lower limit. Clearly the range of scales 
probed is roughly the same as for the large-angle case: Imax/lmin ^ N. Note 
however that one is not assured that one has the complete picture for these 
scales since one also is limited in temporal scales as well (see fig ||). One could, 
in principle, piece together all scales by combining simulations with different 
box sizes and different temporal ranges. Being able to do so presupposes 
that the different angular and temporal scales are essentially uncorrelated. 
This is liable to be less true in non-Gaussian models. For example in string 
models, even long strings at late times which subtend a very large angle, do 
also exhibit temperature discontinuities which will have an effect on all angular 
scales. However the discontinuities from these longer strings are rare and do 
not contribute greatly to the smaller scale anisotropies when compared to the 
much more numerous smaller angle strings. The assumption of independence 
of scales is a good approximation but not exactly correct. Note also that the 
finite box size not only limits the access to large spatial scales but also to 
large temporal scales. Once the horizon sufficiently exceeds the box size, the 
seeds will behave differently than they would in an infinite universe. Thus we 
must also presuppose that there are no long term temporal correlations which 
will be different in periodized box than in an infinite box. Generally we find 
that seeds only produce significant perturbations on scales comparable to the 
horizon size and thus any long term temporal correlations which might exist 
cannot not have much effect on the anisotropies. 

3 The Brightness Distribution 

Here we would first like to point out just how much information is required 
to describe the photon distribution present in a simulation. Recall that in 
a Gaussian models the statistical distribution of the anisotropies observed by 
any given observer is determined completely by the ensemble-averaged angular 
power spectrum, i.e. the Ci's. For angular resolution characterized by /max 
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one only needs ^max real numbers. Note that the C/'s do not give the spatial 
correlations in the anisotropy pattern but, as this is not observable in practice, 
these spatial correlations are not of interest. For a non-Gaussian model where 
one relies on numerical experiment one really wants to use the distribution 
function of anisotropy patterns for all observers in the simulation. If one has 
spatial resolution elements and (/max + 1)^ angular resolution elements 
then one needs naively N^{ltna.x + 1)^ real numbers to describe the brightness 
pattern.^ So, for example, we could describe the brightness pattern as 

/max / 

£ J2 a<,,„,(x)y<,,„,(n;z) . (2) 

1=0 m=-'l 

where the a(i,„)'s are the N^ilmnK + 1)^ real numbers (recall that = 
(— l)™a(;_„)* and there are thus 2/ + 1 real numbers for each I). Since typ- 
ically Zniax ~ N this is numbers - which is a rather formidable amount 
of data. Note that there is some information about the uninteresting spatial 
correlations in these numbers. In fact one can use the spatial correlations to 
significantly reduce the amount of information needed to describe the bright- 
ness distribution. 

The Fourier transform of eq (^) is 

imax I 

EE E «(.™,(k)r<,„,(n;z)e*'' . (3) 

k (=0 m=-i 

which still contains iV'^(/rnax + 1)^ real numbers. Here n gives the direction in 
which the brightness is measured. The 2nd argument of y|i.,„), z, indicates the 
North pole of the spherical polar coordinates used in the spherical harmonic 
decomposition. In this case the m decomposition has no particular physical 
significance. If however we rotate the spherical polar coordinates, differently 
for each k, such that the North pole is the direction of k, we obtain a more 
interesting decomposition: 

AJ^ 'max I 

— (n,x) = ^^ ^ A,,„,(k)r,„„,(n;k)e*^ (4) 

k /=0 m=-l 

"We use (/max -I- 1)^ for the number of of angular resolution elements which is really only 
appropriate for a large-angle simulation, although the number of independent lines-of-sight 
is ~ N'^ for both small and large-angle simulations - and the size of the computation is the 
same in both cases. 
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Mathematically this rotation is given by 



(k)= ^ i^5„^„,(z,k)A,„„,(k) , 



(5) 



where the D', 



m.m' 



I are Wigner rotation matrices defined by 



/ 



d2nr(*,„,(n;z)r(,,„,,(n;k) 



(6) 



which obey the sum rule 



Dmmi (^i Dmm2 (^J ^) ^ 



'mim2 ■ 



(7) 



m— — l 



The m decomposition of eq ^ is useful is because the different m's correspond 
to the different helicities of the brightness pattern. In linear theory (i.e. small 
perturbations) the different helicities of the brightness pattern only couple to 
different components of the gravitational field. In particular m = couples to 
scalar perturbations, m = ±1 to vector perturbations, and m = ±2 to tensor 
perturbations. Helicities with |m| > 2 are not coupled to the gravitational field 
at all. Any brightness anisotropy with |m| > 2 helicities, with no gravitational 
"forcing terms" , would be damped to zero by scattering of the photons. Thus 
we expect that the brightness pattern in our universe to only have helicities 
with \m\ < 2, i.e. scalar, vector, and tensor modes. It would indeed be inter- 
esting, as a check of this notion, to decompose the brightness pattern in our 
patch of the universe into helicities. However this requires a knowledge of the 
spatial pattern of anisotropics, and we cannot expect to measure this anytime 
soon. Nevertheless, in a simulation we do know the spatial pattern, and we 
may use the helicity decomposition, as a form of data compression by ignoring 
the A(, ,„)'s with |m| > 2 which we know a priori to be zero. By doing so we 
have reduced the amount of data required to describe the brightness pattern 
to blmaxN^ ~ N'^, i-e. reducing the dimensionality by one. The special "helic- 
ity state" of the CMBR brightness pattern has long been understooctd and is 
implicit in most anisotropy computations, but the transformation between he- 
licities and the spatial brightness distribution is rarely used as it is not needed 
in Gaussian theories. 

In Gaussian theories the brightness pattern is described by the C;'s and 
in non-Gaussian theories the C/'s still exist but do not give the entire picture. 
Since the C/'s give the expectation value of the |a(, averaged over all real- 
izations one cannot determine their values from a single simulation. However, 
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just as we do on Earth, each one can construct an unbiased estimator of the 
Ci for each observer position in the simulation 



^Kx) = E ^1^. (8) 

m=—l 

and, assuming ergodicity, one should take a volume average of the estimators 
from all the observers: 

k m=— Z 

This is the best estimator one can obtain from a single simulation. To obtain 
the 2nd equality in this equation the orthonormality of the e^^'^ and the sum 
rule for D^^^i^s has been used. Since only scalar, vector, and tensor modes 
contribute one may rewrite this as 



di = cf + cr + c^ (10) 



where 



^ ^ |A(,,o)(k)| 

' 07 4- 1 



k 2^+^ 

A,,+„(k)|2 + |A,,_„(k)|2 



k ^' + ^ 

^ _ ^ |A,,^.)(k)p + |A,,_,,(k)p 

k 

Note that there are no cross-terms from the different helicities, just as there is 
none for different k- modes, which results from the fact that these modes are 
orthogonal. Note that if one only knows the pattern of anisotropy as viewed 
from a given point one does not have enough information to decompose the 
pattern one observes into it's scalar, vector, and tensor parts. One can do so 
only if one knows the pattern throughout space, which we do in the case of a 
simulation. 



4 The Defects 

Usually at late enough epochs and large enough scales to be of interest for 
CMBR anisotropy that the defects only interact with the rest of the mat- 
ter gravitationally. For example in the case of (non-superconducting) cosmic 
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strings the GUT scale strings are very thin and the "scattering width" ex- 
tremely small by astronomical standards. Since the gravitational field of an 
object is completely determined by it's stress-energy tensor it is suSicient to 
describe the evolution of defects by the time-history of their stress-energy ten- 
sor, i.e. 0^^(x, t). This real symmetric matrix has ten components, however 
in order for it to be compatible with General Relativity it must obey energy 
and momentum conservation which reduces the number of independent com- 
ponents to six. 

Defects are usually treated in the "stiff source" approximation which ba- 
sically means that the evolution of the defects are evolved in the metric of 
the unperturbed universe. The reasoning is that since any workable model 
of defects can only produce small metric perturbations which implies that the 
stress-energy of the defects must be small compared to that of the other matter 
in the universe. Taking into account the "back-reaction" of this small pertur- 
bations on the stress-energy of the defect would be to include terms of 2nd 
order in smallness. In mathematical language this means that one expands 
everything to 1st order in 0,,,/. Strictly speaking it is true that some defects 
are highly localized in space and can locally dominated the stress-energy, but 
such defects also have enormous internal stresses, comparable to their energy 
density, and their motion is little effected by the weak gravitational fields of the 
inhomogeneities over the course of a dynamical time. The stiff-source approx- 
imation does miss out on the cumulative effects of back-reaction over many 
dynamical times, such as the decay of cosmic strings into gravity waves, but 
these effects are usually only important on small scales and do not have much 
direct effect on CMBR anisotropics. 

The stiff soiircc approximation combined with linear perturbation theory 
means that one can evolve the defects independently of the rest of the matter 
and that the everything besides this evolution is linear. The equations for 
the evolution of the background matter are linear in Q^i, so the solution will 
consist of a term which is linear in Q^,^ and one which depends on the initial 
conditions. For example the GMBR anisotropy may be written as 



— (n,x,r,)=/ rfr?' / dVZ)''-(n,x-x',r?,77')0M-^(x',r;). (13) 



Here we have used conformal time drj = dt/a where a{t) is the cosmological 
scale factor. The 1?^'^ are given by solutions of the coupled Einstein-matter 



AT _ AT^ AT 

rj-i rj~i I rj-i 



I 



(12) 



where the 1st term is linear in O 
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equations of motion and will depend on the matter content of the universe but 
not at all on the defects. Causality assures us that Z?^" — when |x — x'| > 
7] — rj' . The (AT/T)^ term is often referred to as the compensation. We should 
stress that this S/I decomposition as well as the exact form of D^"^ will depend 
on how one chooses to solve the equations. In particular by substituting any of 
the equations of energy-momentum conservation into eq ([l3| ) and integrating 
by parts one would obtain a different expression for D'"^ and one might also 
incorporate the resulting boundary term on the initial hypersurface into the 
compensation, obtaining a different expression for (AT/T)^. There is, however, 
some physics in the compensation; namely the assumption that the universe 
was completely homogeneous before the production of the defects; at least to 
the extent that the inhomogeneities produced subsequently by the creation 
and motion of the defects dominates over any initial inhomogeneity. We do 
not discuss compensation further as this is not a large part of the computation. 

4- 1 Scalar- Vector- Tensor Decomposition 

We have decomposed the brightness pattern into scalar, vector, and tensor 
parts, and one can do the same for both the stress-energy tensor and the 
Green function D^" . For each k-mode define a orthonormal basis 

I - k 

k m! rhr where k = — . (14) 

The orthonormality determines m[. and rhj^ up to a rotation about k, and we 
do not to specify it any further than that. We may decompose the stress-energy 
tensor as 

Soi = +e"k + eimi + e"m- 
e,, = le°i + e""(kk-ii) 

-He"l(kml +mlk) + e"~(km^ +m~k) 

-|-0^(m'm' — m^m^) + 6^ (rh'rh^ + rh^m') , (15) 

where we have used Latin indices, i,j,k,... to denote spatial coordinat;es. 
Symmetry arguments tell us that the Fourier transform (wrt x — x') of D^^'^ 
may be decomposed as 

i3°'(n,k,77,r;') = £)" k + 1?-^ (cos (/)kml + sin (/)km") 
5*^ (n,k,?7,r;') = + (ick- il) 

4-_D"^ ((km' -t- m'k) cos 0k + (m^k -|- krh~) sin (/)k) 
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^(m'rh'— m m ) cos 2(/)i(; + (m'm +m m')sin2(/)k 

(16) 

where we use the notation 

T r ' f ± " "ii! . n rhr 

I^dij /^k = n-k cos 0k = , ;) sm 0k = , .> ■ (17) 

V 1 - Aik Vl - Mk 

Here the arguments of D"^ arc (/ik, Ikj,??,??'), and there is no dependence on 
</>k- 

4. 2 Anisotropy From the Defects 

Combining eqs ([T^fc|l^) we find that eq (|l^) becomes 

+ +25"-^e"l) cos0k+ (5^6^ +2I)"^e"") sin 0k 



-2D-L-Le+ cos20k + 215^-^6'' sin20k 



(18) 



The above expressions for (AT/T)^ is written so that the 1st, 2nd, and 3rd 
lines contain only scalar, vector terms, and tensor terms, respectively. Since 
this equation represents a purely gravitational effect it contains no higher he- 
licities. If one compares eq ( p^ ) to eq (^ and defines cos6'k = /^k one sees 
that the angles (0k,0k) are just the arguments of the y(, ,„)'s in the helicity 
decomposition of the brightness pattern. By equating terms one see that 

A,,o,(k) = z'\/2^ dri (dt^,, + ii^J'e" + + ^i?l'"e"" 
A,,±„(k) = -^'^f H (©' T *e-) + 2D\^ (e"i t *e" 

A(,.±3)(k) =^'V2^ / d77'A^^(e+T«e^) (19) 
where we have expanded the -D's in terms of spherical harmonics: 

_ oo ^ { ^ ^=00,11,0,1111 

i?^ = [7;](-l)'"y2^^z'Y;,,„,(0k,O)A-^ m=\\ x=±,i,± (20) 



2 X-- 
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The extra factor, [i], is only for i^H and D-^ since with this definition all the 
Df's are real. 

We see from the above that there are at most 7 Green functions which 
one must compute, 4 scalar {D°°, D\ D°, D"), 2 vector {D-^, i)"-^), and 
one tensor (D^^). However one can use the equations of energy- momentum 
conservation, which may be written 

§00 + -(Ooo + e°) = ike' + 2-e" = ik{@° + ^e") 

a . a 3 

e'- + 2-e- = ifceii- ei + 2-ei = ifce"i (21) 

a a 

to eliminate 2 of the scalar terms, and 1 (for each of | and — ) of the vector 
terms. Thus one really only needs to calculate 4 Green functions, 2 scalar, 
1 vector, and 1 tensor. The vector and tensor Green functions are each used 
twice, once for each "polarization" of vector and tensor modes. 



4-3 Computational Cost 

Now consider the cost to compute the full anisotropy pattern from a simulation 

of the defects. We do not include here the cost of simulating the defects, i.e. 
producing the Qni,{x, rj), which typically scales as iV. For a cubical simulation 
box the Fourier spectrum take on the discrete values 

k = (22) 

where L is the size of the cube, and the components of the vector N is any 
triplet of integers. With spatial resolution elements there arc N"^ modes 
and the components of N take on the values {—N + 1, —N + 2, . . . , iV — 1, N}. 
The Green functions depend only on |k| which is 27r/L times the sum of the 
squares of 3 integers in this range. Clearly there are no more than 3(7V/2)^ 
such values. For a million different k's there are only a few thousand different 
|k|'s. Usually one would want to compute these Green functions for one fixed 
observation tim(\ r/, but would need to have temporal resolution in rj' the same 
as in the simulation. The Courant condition suggest that this resolution must 
correspond to at least timesteps during the simulation. Since the Green 
functions are usually small ii I ^ k{'q — rj' ) we may a priori ignore terms with 
larger I. Since the vast majority of modes have k ^ N/ and most of the 
required terms have (ry — rj') > L; we see that typically we want Di's for I 
up to /max ^ N. Thus a complete table of Green functions would have size 
N^. One could certainly reduce this greatly since, in principle, one needs 
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a wavenumber resolution no better than Afc ~ 27r/L, reducing the number 
of fc-modes to something of order N rather than iV^. To compute each of 
these Green functions numericaUy requires in integration from rj' to ry, taking 
roughly N operations; so the cost of computing these Green functions scales 
like A''"' if one does not oversample in k {N^ if one does). 

Performing the integrals of eq (|l^) is liable to be a larger task, as this 
clearly scales as N^; there being different k's and 5?max integrals for every 
k and the integrals themselves requiring ~ N operations. We do stress that 
the prefactors multiplying these scaling laws can be reduced significantly by 
not computing every I or coarsening the resolution of the integrations. Once 
one has done the integrals one has the A(,^)'s from which one could compute 
the best estimates of the angular power spectrum using eq (|^) , which is itself is 
an N'^ computation. However the power spectrum is not the entire story and 
one would really like to have the full anisotropy pattern for of the observers in 
the simulation volume. To do this one needs to perform the Wigner rotation 
of eq (||), which converts the ~ iV* A(, ,„)'s to the ~ a(i,m)'s and is an 
A'^^ computation. Finally to transform the 's to a(, ,„)'s one may use an 
FFT which is again an computation. Note that if one does want to know 
the full brightness distribution then the final data set is ^ in size and the 
computational cost scales in the same way. Clearly it could not scale in any 
better way. 



5 Computational Reduction 

We have seen that the cost of computing the full brightness distribution from a 
defect simulation scales as N^, which is more costly than the cost of simulating 
the defects (iV'*), and in any case makes the computation large. How can one 
can get away with doing less? For computations of large-angle anisotropy the 
effects of scattering is not important and the anisotropy in defect models is 
largely determined by the Sachs- Wolfe integral along the line-of-sight, perhaps 
with the inclusion of a boundary term at the last-scattering surface. In this 
case the brightness along a given line-of-sight is independent of the brightness 
along other lines-of-sight and one may compute these separately if one knows 
the metric perturbation throughout space-time (an TV* computation). Thus 
one can con* 
of observers! 



S effort on the photons converging on a small number 
performing an computation for each observer .| 
When scattering is important, as for small- angle anisotropy, the brightness 

''Bouchet et al. i and Allen et al. use a somewhat different technique which avoids com- 
puting the metric perturbation directly, and is particularly advantageous for cosmic strings 
which do not fill space. 
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along a given linc-of-sight will, strictly speaking, depend on the brightness in all 
other directions at all points along that lines-of-sight. Since the number of these 
brightnesses is ~ N^, one might worry that one requires an computation 
to find the anisotropy pattern for just one observer! However in practice one 
finds that this dependence is almost exclusively on the first few moments of the 
brightness pattern (e.g. monopole and dipole) and depends only very weakly on 
the higher moments, and that one can compute these lower moments without 
much regard to thcJiigher moments. These sorts of approxinmtiDns have been 
in use for decadefru they have been perfected only recentljE^lo. In practice 
computing the anisotropy with scattering is no larger a computation than 
without, one need only do a few integrals along the line of sight! Note however 
that if one were to attempt computing the anisotropy pattern in this way 
for all N'^ observers in the box one would end up with a computation which 
scales as ~ N^. Given that the results for different observers from the same 
simulation are not completely statistically independent it is not necessarily 
worthwhile to compute the anisotropy— pattern for all observers allowed by 
ones spatial resolution. Allen et al. Illl have used 64 observers in a single 
large-angle simulation and found no significant correlations between observers. 
Given that the cost per observer in a given simulation goes as N'^ while the 
cost of a simulation goes as it is certainly worthwhile to use many observers 
in a given simulation. 



6 Another Algorithm 

When there is scattering one may restrict oneself to a few ohaepzErs but still 
use the Fourier representation described above. The classic us(alZlli3 of a small- 
angle simulation is to only compute the brightness for lines-of-sight traveling 
along one or more of the 6 principle axes of the box (±x, ±y, ±z), as illustrated 
in fig |2. To do this one need just use the appropriate values of /ik and (j)]^ 
in eq (18|) for the direction one is interested in. What one obtains is the 
brightness pattern, AT/T{h, x), but with fixed h. As long as h is along one of 
the principal axes then each slice of this function of x perpendicular to fi gives 
the brightness pattern which would be seen by a "distant" observer on a square 
patch of sky. There would be N different slices which represent the sequence 
of patterns a given observer would see over a period of time. Performing the 
integrals is an ~ N'^ algorithm (FFTing to real space only ~ N'^) so the cost 
per pattern is ^ N'^, the same as required for a similar result using the line- 
of-sight method mentioned above. Of course there will be some correlations 
between the patterns on nearby slices. 
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7 Summary 



This report lists some of the mathematical details of CMBR anisotropy which 
could be used in computations of CMBR anisotropy from defects. Hopefully 
these details are interesting even to those not planning to do the calculations. 
Many other details were not mentioned for lack of space, specifically those 
relating to computing the Green functions. The equations used are just the 
usual cosmological Boltzmann-Einstein equations with the appropriate initial 
conditions. The formalism specifically outlined has not been used in the past 
(which is why it was concentrated on) but some papers |applying some of these 
technique to cosmic string simulations will appear sooirJ. 
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